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We report a comprehensive micro-Raman scattering study of electrochemically-gated graphene 
field-effect transistors. The geometrical capacitance of the electrochemical top-gates is accurately 
determined from dual-gated Raman measurements, allowing a quantitative analysis of the frequency, 
linewidth and integrated intensity of the main Raman features of graphene. The anomalous behavior 
observed for the G-mode phonon is in very good agreement with theoretical predictions and provides 
a measurement of the electron-phonon coupling constant for zone-center (T point) optical phonons. 
In addition, the decrease of the integrated intensity of the 2D-mode feature with increasing doping, 
makes it possible to determine the electron-phonon coupling constant for near zone-edge (K and 
K’ points) optical phonons. We find that the electron-phonon coupling strength at T is five times 
weaker than at K (K’), in very good agreement with a direct measurement of the ratio of the 
integrated intensities of the resonant intra- (2D’) and inter-valley (2D) Raman features. We also 
show that electrochemical reactions, occurring at large gate biases, can be harnessed to efficiently 
create defects in graphene, with concentrations up to approximately 1.4 x 10 12 cm -2 . At such 
defect concentrations, we estimate that the electron-defect scattering rate remains much smaller 
than the electron-phonon scattering rate. The evolution of the G- and 2D-mode features upon 
doping remain unaffected by the presence of defects and the doping dependence of the D mode 
closely follows that of its two-phonon (2D mode) overtone. Finally, the linewidth and frequency of 
the G-mode phonon as well as the frequencies of the G- and 2D-mode phonons in doped graphene 
follow sample-independent correlations that can be utilized for accurate estimations of the charge 
carrier density. 

PACS numbers: 78.67.Wj, 78.30.-j, 72.80.Vp, 63.22.Rc, 63.20.kd, 82.45.-h 


I. INTRODUCTION 

Graphene, as an atomically thin two-dimensional crys¬ 
tal, features an electron gas that is directly exposed to 
its local environment. As a result, graphene is uniquely 
sensitive to external stimuli. This is remarkably illus¬ 
trated by the electric field effect, which makes it possible 
to swiftly tune the carrier density of graphene (i.e., its 
Fermi energy Ep) and, in return, to control a wealth of 
fundamental properties, among which, the electrical 1 3 
and optical conductivities,®® as we ll as of the electron- 
phonon coupling? 8 From a more applied standpoint, the 
unique controllability of graphene can be harnessed in a 
variety of nano-devices.® 

Among the various experimental techniques empl oyed 
to study graphene, Raman scattering spectroscopy® 11 
stands out as a fast, sensitive, and minimally invasive 
tool in order to probe electron-phononf 7 * 8 * 121 ®! electron- 
electron 15 and electron-defect scattering 16 ^ at vari¬ 
able carrier density. Raman spectroscopy is also rou¬ 
tinely e mploye d to characterize unintentional doping in 
graphene®*® an d to study the sensitivity of graphene 
to atmospheric 21 and chemical dopants. 22 28 Quantita¬ 
tive investigations of doped graphene are particularly 
relevant, since several interesting phenomena, such as 
superconducti vity)® !® ferromagnetism^ 2 charge or spin 
density waves ,® 1 ® as well as changes in the plasmon 
spectrunf®® are expected to occur in the strong dop¬ 


ing regime (|-Ep| > 1 eV). 

In practice, solid state graphene field-effect transistors 
(FETs), typically using a Si substrate as a back-gate and 
a Si 02 epilayer as a gate dielectric, have been widely 
used to study the Raman response of graphen e in the 
vicinity of the Dirac point (|-Ef| = 0 — 300 meV)! 7 ® 7 ® 
To access higher dopin g l evels, other methodologies 
based on chemical doping 22 "® an d electrochemical gat¬ 
ing 39 41 have been introduced. The former is highly 
efficient, resulting in charge carrier concentrations ex¬ 
ceeding 10 14 cm -2 , but is irreversible and little con¬ 
trollable. The latter, which relies on the forma¬ 
tion of nanometer-thin electrical double layers (EDL) 
with high geometrical capacitance, makes it possible 
to reversibly attain electron or hole concentrations 
as high as ~ 10 14 cm -2 , (i.e., |JTp| ~ 1 eV) at cryo¬ 
genic temperatures!® Recently, electrochemically-gated 
graphene FETs have been succe ssfully employed to inves¬ 
tigate electron-phonon coupli ng! 41 1 4 3 " ® but also bandgap 
formation in bilayer graphe ne)®- 49 electron transport at 
high carrier density)®®® many-body phenomena,® as 
well as to electrically contr ol th e interaction between 
nano-emitters and graphene! 52 * 53 ! 

In such studies, an accurate determination of Ep 
(hence of the gate capacitance) as a function of the gate 
voltage is a critical requirement. However, as opposed 
to solid state FETs, in which the oxide dielectric con¬ 
stant and thickness can be known with accuracy, the 
thickness of the electrical double layer may be highly 
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sensitive to the device geometry, fluctuate spatially and 
vary over time. This further highlights the need for 
(i) robust methods for device fabrication and, (ii) ac¬ 
curate tools to experimentally measure the gate capaci¬ 
tance. In previous works, the gate capacitance of elec- 
trochemically gated graphene FETs has been evalu ated 
from an estimation of the thickness of t he EDL ,* 41 * 44 * from 
capacitanc e! 54 * 55 ! or Hall measur emen ts, ^ 42 * 50 * 51 * or from 
optical absorption spectroscopy.*^*^* 

In this article, we show that micro-Raman scattering 
measurements on electrochemically top-gated and SiC >2 
back-gated graphene FETs can be used to accurately de¬ 
termine the geometrical capacitance of the electrical dou¬ 
ble layer and hence Ep, with a spatial resolution down 
to approximately 1 /iin. Calibrated electrochemical gates 
allow us (i) to quantitatively compare the anomalous 
doping-dep enden ce of the G-mode phonon to theoreti¬ 
cal modelsj^^’i! (ii) to deduce the electron-phonon cou¬ 
pling constants at the center (T point) and near the ed ges 
(K and K' points) of the Brillouin zone of graphene*^* 
and (iii) to establish well-defined correlations between the 
frequencies, linewidths and integrated intensities of the 
main Raman features in doped graphene. Importantly, 
we show that at top-gate voltages beyond the threshold 
for electrochemical reactions, defect concentrations of up 
to approximately 1.4 x 10 12 cm -2 can be created with¬ 
out damaging the device. This allows us, in particular, to 
quantitatively investigate the doping dependence of the 
defect-related D mode and to estimate the electron-defect 
scattering rate in graphene. 

The paper is organized as follows: the experimental 
methods are exposed in Sec. [TTJ Section [TTT] presents a 
model for the electric field effect in electrochemically- 
gated graphene FETs. Section [TV] is dedicated to the ex¬ 
perimental determination of the geometrical capacitance 
of the electrical double layer. In Sec. [V] we specifically 
address electron-phonon coupling in pristine graphene. 
Section |Vl] describes our study of defective graphene and 
the determination of the electron-defect scattering rate. 
Finally, in Sec. |VII| we describe the correlations between 
the main Raman features in doped graphene. 


II. METHODS 
A. Sample preparation 

Graphene samples are produced by mechanical exfolia¬ 
tion of natural graphite onto highly p-doped Si substrates 
covered with a (285 ± 15) nm SiC >2 epilayer. Graphene 
monolayers are identified by optical microscopy and 
micro-Raman spectroscopy. Source, drain and gate elec¬ 
trodes are made by photolithography, followed by metal 
deposition (Ti (3 nm)/Au (47 nm)). The device are then 
coated with a ~ 4 (im thick photoresist layer (MicroChem 
SU8 2005), and a second photolithography step is per¬ 
formed to open a window above the graphene channel 
and gate electrode, as shown in Fig.[l](a)-(b). Finally, the 


electrochemical top-gate is formed by depositing a drop 
of polymer electrolyte with a micropipette. The poly¬ 
mer electrolyte is prepared by mixing lithium perchlo¬ 
rate (LiCICH) an d polyet hylene oxide (PEO) in methanol 
at a weight ratio 17 ^^* 0.012:1:4. The mixture is then 
heated at 45 °C and stirred until it becomes uniform. 
This suspension is filtered to get a clear solution. Af¬ 
ter dropcasting, the methanol evaporates and a thin film 
of transparent polymer electrolyte is formed. To remove 
residual moisture and solvent, the devices are annealed at 
about 90 °C. Noteworthy, the device geometry depicted 
in Fig. [l](a)-(b) features a well-defined gated region and 
prevents the polymer electrolyte to be in contact with 
the so urce an d drain electrodes. As compared to earlier 
works?- 41 * 43 * 44 * this improves the gating efficiency and re¬ 
duces the electrochemical reactivity of our devices. Some 
measurements described in the following are performed 
in a dual-gated geometry. In this case, the back-gate 
voltage is applied using the Si substrate as a gate elec¬ 
trode. 


B. Experimental setup 

We perform micro-Raman scattering measurements 
in ambient conditions on top-gated and dual-gated 
graphene field-effect transistors. Raman spectra are 
recorded in a backscattering geometry, with a home-built 
setup, using a x 40 objective (NA = 0.60) and a 532 nm 
laser beam focused onto a spot of approximately ~ 1 /irn 
in diameter. The sample holder is mounted onto a x- 
y-z piezoelectric stage, allowing spatially resolved Ra¬ 
man studies. The collected Raman scattered light is dis¬ 
persed onto a charged-coupled device (CCD) array by a 
single-grating monochromator, with a spectral resolution 
of about 1 cm -1 . The laser beam is linearly polarized and 
the laser power is maintained below 500 /iW, in order 
to avoid thermally induced spectral shifts or lineshape 
modifications of the Raman feat ures,^ as well as photo¬ 
electrochemical reactions? 16 * 42 * 45 * The sample holder is 
electrically connected to a sourcemeter, which triggers 
our CCD array. Raman spectra are recorded as a func¬ 
tion of the applied gate bias, once a steady gate leak 
current (typically lower than 100 pA in the electrochem¬ 
ically top-gated configuration) is achieved. For this pur¬ 
pose, the gate bias is first applied for a settling time of 
~ 1 min, before recording each Raman spectrum. This 
procedure ensures that Raman spectra are recorded at 
constant charge carrier densities. Raman spectra are also 
recorded during several forward and backward top-gate 
sweeps at the same spot on a given sample and very 
reproducible results, with no significant hysteresis, are 
observed. We find, however, that the geometrical ca¬ 
pacitance of the top-gate, as well as the electron-phonon 
coupling constant may exhibit a certain degree of spatial 
inhonrogeneity. Additionally, in ambient air, the gate ca¬ 
pacitance may decrease over time, by up to one order 
of magnitude over a couple of days, due to a degrada- 
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(a) (c) 



Polymer electrolyte ^—Photoresist ^—Graphene 


FIG. 1. (a) Optical image of a dual-gated graphene field-effect transistor prior deposition of the polymer electrolyte. The 
source and drain electrodes are covered with photoresist (SU8) to prevent them to be in contact with the polymer electrolyte, 
(b) A schematic cross-section of our dual gated graphene field-effect transistor, with Li + (green) and CIOJ (red) ions and the 
electrical double layers near each electrode. The Si substrate is used as a back-gate, (c) Schematic energy diagrams of the 
electronic states of the gate electrode and of graphene. Occupied states are represented in grey. At zero gate voltage (V = 0), 
the electrochemical potentials of the gate electrode p and the graphene layer /.aslg are equal. The Fermi energy of graphene is 
Ep. Applying a gate voltage V results in an electrostatically-induced shift ( ecj )) and in a change of the Fermi energy of graphene 
(Ef). The electrochemical potential difference is equal to eV, leading to Eq. ([!]) with eVo = Ep. (d) Equivalent electrical circuit 
of our device at steady state. Cbg is the geometrical capacitance of the Si/Si02 back-gate, Cg is the geometrical capacitance 
of the electrical double layer at the graphene/polymer electrolyte interface and Cq is the quantum capacitance of graphene. 


tion of the polymer electrolyte. Such aging effects un¬ 
derscore the necessity of fast characterizations of electro- 
chemically gated FETs and may account for the fairly 
large spread in the gate capacitances reported in litera¬ 
ture. In order to avoid sample aging effects, our mea¬ 
surements were performed immediately after deposition 
of the polymer electrolyte. Interestingly, the dispersions 
obtained from a set of measurements at several spots on 
a given graphene FET are very similar to the sample- 
to-sample dispersions observed by measuring at (single) 
random spots on a set of graphene FETs. This further 
highlights the interest of spatially resolved studies. 


III. ELECTRIC FIELD EFFECT 

Figure [2] shows typical Raman spectra recorded over a 
top-gate voltage sweep, with the two prominent Raman 
features in pristine graphene: the first order G-mode fea¬ 
ture, which involves zone-center optical phonons (at the 
r point), and the second-order resonant 2D-mode fea¬ 
ture, which involves near z one-edge optical phonons (at 
the K and K’ points).^ ^ Note that no defect-induced 


D-mode feature emerges from the background in our ex¬ 
perimental conditions. This illustrates the high struc¬ 
tural quality of the graphene sample. As expected,^ the 
G-mode frequency and linewidth vary significantly with 
the top-gate bias (Vtg)- Similar trends are observed by 
applying a back-gate voltage (Vbg)- The minimum value 
of the G-mode frequency wq and the maximum value of 
its full width at half maximum (FWHM) Tg are reached 
at the same value of Vtg.o = +0.5 V. This value cor¬ 
responds to the charge neutrality point (CNP), where 
Ep = 0. The CNP is reached at a finite Vtg.Oi due to an 
unintentional doping of the graphene layer, induced by 
the substrate as well as the polymer electrolyte.^ A fi¬ 
nite value of Vtg — Vtg,o results in a finite charge carrier 
density n. In this work, a positive (negative) gate voltage 
corresponds to electron (hole) injection.^ Qualitatively, 
for both positive and negative values of Vtg — Vtg.ch we 
observe a nearly symmetric increase of wg accompanied 
by a symmetric decrease of Tg (see Sec. V A for details). 
In contrast, the 2D-mode feature is less sensit ive to dop¬ 
ing than the G-mode featur^B (see Sec. VB). 

In order to carefully study the G- and 2D-mode fea¬ 
tures as a function of Ep, one has to convert the gate 
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FIG. 2. (a)-(c) Color maps of the Raman spectra of a pristine graphene monolayer (sample 1), measured using a 532 nrn laser 
beam, as a function of the top-gate voltage Vtg- The G- and 2D-mode features appear prominently and no defect-induced 
D-mode feature is observed. Panels (b) and (c) show a clear evolution of the G- and 2D-mode features with varying Vtg. 
The black dashed lines correspond to the central frequency of each Raman feature. The charge neutrality point (CNP) is 
indicated by an arrow, (d) Ranran spectra at values of Vtg between -0.5 V and +1 V. The circles are the experimental data 
and the solid lines are fits (see text for details). The CNP is reached at Vtg.o = +0.5 V (see green line), (e)-(g) Schematic 
representation, in the momentum-energy space, of the main Raman features in graphene! 1 TI (hus) denote the incoming 
laser (Raman scattered) photon energy. The G mode (e) is a non-resonant process.^ 5 46 In (f), the 2D mode is represented 
as the dominant fully resonant inner process p ,6t) involvi ng two near zone-edge transverse optical phonons with frequency cud 
and opposite momentum ±q. The defect-induced D mode 11 60, «J (g) involves one electron-phonon (solid arrow) scattering (q, 
cjd) and one electron-defect (dashed arrow) scattering process. One of these two processes is resonant. In (g), a resonant 
electron-phonon scattering process is represented. 


voltage into Ep or, equivalently, n. First, the Fermi en¬ 
ergy at a given n is Ep = sgn (n)HvF^n\n\, where H is 
the reduced Planck’s constant and vp ss 1.1 x 10 6 m/s 
is the Fermi velocity of graphene on a SiC >2 substrate.^ 
Note that this formula applies only at T = 0. How¬ 
ever, in practice, finite temperature effects only induce 
a very minor correction to this simple scaling!^ 5 An ap¬ 
plied top- or back-gate voltage V creates an electrostatic 
potential difference <j> between the graphene monolayer 
and the gate electrode. Besides, the injection of charge 
carriers in graphene leads to a shift of its Fermi energy. 
Consequently, V introduces a difference in the electro¬ 
chemical potentials of the gate electrode fi and of the 


graphene layer /z S lg (see Fig.JTJc)) 

Mslg — A* = eV = E f + ecf> + eV 0 , (1) 


where e is the elementary charge, Vo is a constant that 
accounts for the initial doping and implicitly includes the 
work function difference between the two materials.^ 
Assuming that the gate can be modeled as a parallel 
plate capacitor with a geometrical gate capacitance Cq, 
the relation between V and Ep is given by 


v_v 0= A +sgn(£F) _^a_. 


( 2 ) 
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Importantly, the first term on the right hand side of 
Eq- © scales as Ep (i.e., ^/n) and is related to the quan¬ 
tum capacitance of graphene Cq P while the second term 
scales as Ep (i.e., as n) and is related to the geometrical 
gate capacitance C G l 41 l 44 l 

For a typical SiC >2 back-gate insulator, the geometrical 
capacitance Cbg per unit area is simply given by Cbg = 
£ r Eo/d,BG > where e r ~ 4 is the relative permittivity of 
SiC> 2 , Eq the vacuum permittivity and c?bg is the Si02 
thickness. In this work, dpG = (285 ± 15) nrn results in a 
back-gate capacitance Cbg = (12.4 ± 0.7) nF cm -2 . For 
a typical Fermi energy Ep ~ 100 meV, the quantity Ep /e 
is negligible as compared to the other term in Eq. §• 

The case of the polymer electrolyte top-gate is slightly 
more complicated. Indeed, when a voltage is applied be¬ 
tween the gate and the SLG, Li+ and CIOJ diffuse in 
the polymer to form electrical double layers at the in¬ 
terfaces as it is sketched in Fig. [TJb). 41 These EDL can 
be modeled as parallel plate capacitors with a thickness 
given by the Debye length dpG, and a geometrical capac¬ 
itance per unit area Ctg = £r£0 /drG ■ The total geomet¬ 
rical capacitance of the polymer electrolyte is thus given 

by Ctg (Sp-gate + ^p-graphene) » wh ere S'p-gate (resp. 
iSp—graphene) is the contact area between the polymer elec¬ 
trolyte and the gate electrode (resp. the graphene mono- 
layer). Since S'p-gate > S p _g ra phene (see Fig. [ija)), one 
only needs to take into account the geometrical capac¬ 
itance of the EDL at the graphene-polymer electrolyte 
interface. The Debye length is theoretically given bjjm 
<7 tg = 2Ce 2 /soErk^T, where T is the temperature, fcs is 
Boltzmann’s constant and C is the concentration of ions 
in the polymer electrolyte. In practice, the exact value 
of C cannot be measured. One can nevertheless obtain 
an estimate of Ctg ~ 4.4 [iF cm -2 , assuming a typical 
value of g?tg ~ 1 nm and £ r ss 5 for PEOpJJ This ca¬ 
pacitance is more than two orders of magnitude larger 
than Cbg and becomes comparable to the quantum ca¬ 
pacitance for Ep ~ 100 meV. As a result, the two terms 
in Eq. § are of the same order of magnitude and must 
be taken into account in the present study. 


IV. GEOMETRICAL CAPACITANCE OF THE 
ELECTRICAL DOUBLE LAYER 


(a) 




(c) 


V TG (V) 



Our first objective is to precisely determine Ctg- Pre¬ 
vious works on oxide dual-gated graphene FET^Ml 
have shown that provided one geometrical capacitance 
is known, the other can be determined by monitoring 
the minimum (source-drain) conductivity point as a func¬ 
tion of the bottom and top-gate biases. At steady state, 
our dual-gated graphene FETs have the same equiva¬ 
lent electrical circuit (see Fig. 0d)) as the devices of 
Ref. [69] Here, rather than using electron transport 
measurements, we apply micro-Raman scattering spec¬ 
troscopy, which provides a local measurement. For a 
fixed VbGj we sweep Vtg and record Raman spectra, 


FIG. 3. (a) Frequency cjg and (b) relative FWHM AFg of the 
G-mode feature as a function of the top-gate voltage, recorded 
at various back-gate voltages on sample 1. The curves are 
vertically offset by 10 cm -1 for clarity. The symbols are 
experimental data, (c) Top-gate voltage Vtg, neutral, corre¬ 
sponding to the CNP in dual-gated graphene, as a function 
of the applied back-gate voltage Vbg- A top gate capacitance 
Ctg = 3.3 fj,F cnC 2 is deduced from a linear fit of the data 
(solid line). The solid lines in (a) and (b) are fits based on 
Eq. i> and |6j), respectively, with Ctg = 3.3 /rF cm 2 . 
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as described in Sec. IIB Then, we extract uq and Tq 
from Lorentzian fits. Figures [3ja)-(b) show these two 
quantities as a function of Vtg f° r five different val¬ 
ues of Vbg- We observe a clear shift of the CNP, at¬ 
tained at Vtg, neutral, with Vbg- In practice, Vtg, neutral 
is extrac ted fr om the Fg(Vtg) curves, which, expectedly 
(see Sec. VA), exhibit a sharper extremum near neutral¬ 
ity than the wg(Vtg) curves. As shown in Fig. [3](c), 
Vtg, neutral varies linearly with Vbg- Indeed, from the 
equivalent circuit in Fig. [ljd), the total charge density 
injected by top- and back-gates leadtP^ 


ne = -Ctg ( Vtg — Vtg.o-— 


— Cbg ( Vbg - Vbg,o-) • (3) 


At the CNP, n = 0 and Ep = 0. Therefore, 


n (x 10 12 cm 2 ) 


-15.3 -9.8 -5.5 -2.4 -0.6 0.0 0.6 2.4 5.5 



VtG, neutral 


Vtg.o 


Cbg 

Ctg 


(Vbg — Vbg.o)- 


( 4 ) 


Using Eq. Q, a linear fit of the data in Fig. 1°) 
yields Cbg/Ctg = (3-8 ± 0.2) x 10 3 . Since 

Cbg = (12.4 ± 0.7) nF cm -2 , we deduce that Ctg = 
(3.3 ± 0.3) /jF cnV 2 , which is of the same order 

of mag nitude as w hat was reported before for similar 
deviees P f 4 41 l 42 l 4 f lZSI \\' e m ay now convert Vtg into Ep. 


FIG. 4. Frequency o>g (red squares, left axis) and relative 
FWHM AFg (blue circles, right axis) of the G-mode fea¬ 
ture, extracted from the measurements in Fig. [2] as a func¬ 
tion of Fermi energy or doping. The corresponding Feynman 
diagrams are shown as insets. The left inset represents the 
renormalization of the G-mode phonon frequency due to in¬ 
teractions with virtual electron-hole pairs. The right inset 
represents lifetime broadening due to the resonant decay of 
a G-mode phonon into an electron-hole pair. The solid blue 
and red lines are fits based on Eqs. (j5j and §, respectively. 
The fitting parameters are Ctg = 3.9 n F cm -2 , Ar = 0.027 
and 5Ep = 35 meV. 


V. ELECTRON-PHONON COUPLING IN 
PRISTINE GRAPHENE 


is given by 


A. Doping-dependence of the G-mode feature 


Considering only lattice expansion, due to the addition 
of charge carriers, one may expect the G-mode frequency 
to increase (decrease) under hole (electron) dopingP^ 
Thus, the peculiar, nearly symmetric behavio rs observed 
here and previously reported by otherVd - 8 ! 41 l 45 l EAZl contrast 
strongly with the trends predicted if one only considers 
lattice expansion effects. This anomalous behavior has 
been originally predicted by AndcP^ and by Lazzeri and 
Mauri 5 ^ as a consequence of the strong coupling between 
zone-center optical phonons and low-energy electronic ex¬ 
citations across the gapless bands of graphene. Related 
effects occur in metallic carbon nanotubesf-^The anoma¬ 
lous doping dependence of th e G-mode can be described 
using the phonon self-energyP^ 7 — — the real part of 
which is equal to wq and the imaginary part to Tg- As a 
result, the evolution of wq and Tq are deeply connected 
(see Fig. [I]). 

The variation of Tg is due to the decay of the G phonon 
into an electron-hole pair (see right inset in Fig. [4]) and 


Ar G = r G — r 0 


Ar o 

= T^g x 


/ - 


flUn 


-Ep 


( 5 ) 


where Wq is the phonon frequency at Ep = 0, f(E) = 
[1 + exp (E/ksT)]- 1 is the Fermi-Dirac distribution at a 
temperature T and Ar is a dimensionless coefficient corre¬ 
sponding to the electron-phonon coupling strength^ (see 
also Sec. VC). r 0 contains all other sources of broadening 


that are independent on the carrier density (anharmonic 
coupling^- 3 disorder, instrument response function). For 
|Ef| > fiwg/2, ATg vanishes due to Pauli blocking. 

The evolution of loq with Ep is the sum of an adiabatic 
contribution loq corresponding to the modification of the 
equilibrium lattice parameter and a non adiabatic one 
Wq A corresponding to the renormalization of the G-mode 
phonon e nergy due to interactions with virtual electron- 
hole pairs' 5 -^® (see left inset in Fig. [4|. At a finite tem¬ 
perature T, the frequency shift is given b; 


Awg = WG — Wg = AWf 


Aw£ A , 


( 6 ) 
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with 


Aw na = iLi {f(E-E F )~f(E)]E^gn(E) ^ 


2 tt hj_ c 


E 2 - M) 2 /4 


( 7 ) 


where f denotes the Cauchy principal value. One should 
note that Arc and Awe are proportional to Ar- 

In the simulations described below, we use the results 
of the calculation by Lazzeri and Mauri to include adia¬ 
batic contribution A cjq (see Eq. (3) in Ref. EZD- Impor¬ 
tantly, for \Ep\ < 1 eV, the adiabatic contribution pro¬ 
vides only a minor correction to the non-adiabatic term 
and does not affect Tq. 

Moreover, to accurately describe the experimental evo¬ 
lution of the G mode, one also has to take int o account 
random spatial fluctuations of the Fermi energyP^IEI '5 
It is reason able to assume that Ep follows a Gaussian dis- 
tributior l 65 l 74 l 75 l around its mean value, with a standard 
deviation SEp. Thereafter, the computed Acjg(Ff) and 
ATg(-E’f) used to fit our data are given by the convolu¬ 
tion of this Gaussian distribution with Eq. ([5| and ([b]) . 

Figure [3] displays the results of simultaneous fits of 
Awg(I / tg) and ATg(Vtg) f° r h ve top-gate sweeps at 
different Vbg- We used wf = 1-1 x 10 6 m s” 1 and the 
values of Ctg, Vtg, neutral and Wq obtained in Sec. |IV| 
Thus the fitting parameters are Ar, SEp and To- The 
experimental data are remarkably well fitted by the the¬ 
oretical model. Interestingly, although the two phonon 
anomalie d 37 * 56 * 57 ! predicted at Ep = ±huiQ by Eq. (|6j) are 
largely smeared out at room temperature, one can still 
notice a hint of their presence in Fig. 1*0.0 and[9](a). 

From these five fits, we get Ar = 0.036 and SEp = 
40 meV. Since SEp « 50meV on bare Si02 without an 
electrochemical top-gatepE 6 we conclude that charge in¬ 
homogeneity doe s no t have a major effect on our analysis. 
DFT calculation^ have predicted Ar = 0.028, which 
is slightly smaller, but consistent with our measurement. 

Another way to further compare the experimental data 
and theory is to set Ctg as adjustable parameter when 
fitting Awg(Ttg) and ATg(Vtg)- This yields Ctg = 
3.9 /j,F cm -2 , Ar = 0.034 and SEp = 35 meV. These 
values are very consistent with the more constrained fits 
discussed above (see Sec. IV). Similar studies were re¬ 


peated on more than five samples, with similar conclu¬ 
sions. This demonstrates that a direct fit of Auig(Vtg) 
and ATg(Vtg) can be used to get an accurate mea¬ 
surement of Ctg, which allows to convert Vtg into Ep 
through Eq. ([2]). This is a much faster approach to deter¬ 
mine Ctg, which does not require a dual-gated device. 
As an example, a fit of the data in Fig. [2] is shown in 
Tig- 0 and shows a very good agreement between exper¬ 
iment and theory. More generally, our fitting procedure 
allows us to estimate Ctg and Ap with relative uncer¬ 
tainties of approximately 20% and 10%, respectively. 

To better understand the importance to fit simultane¬ 
ously Awg(Vtg) and ATg(Vtg), we have fit these quan¬ 
tities separately for the measurements shown in Fig. [3] 
(not shown). From the fit of Au;g(Vtg), we obtain 
Ctg = 2.3 /./F cm" 2 , A r = 0.042 and SEp = 50 meV. 



FIG. 5. Fermi energy Ep as a function of the relative fre¬ 
quency of the G mode Au;g- Measurements on five different 
devices are represented with different symbols. The dashed 
and solid lines correspond to Eq. § and |To|), respectively. 


Except the large value of Ap, the two parameters are rea¬ 
sonable. From the fit of ATg, we obtain Ar = 0.033 
and SEp ss 40 meV and an unrealistically large Ctg ~ 
100 [iF cm” 2 . The latter value suggests that the behavior 
of ATg can be rationalized using solely the quantum ca¬ 
pacitance of graphene. This is understandable, since the 
variations of ATq occur near Ep = 0, where the contri¬ 
bution of the quantum capacitance dominates in Eq. |2| . 
However, the value of ATq near Ep = 0 is directly pro¬ 
portional to Ar and is not influenced by C TG , while Accq 
varies mostly away from the CNP. Hence, its evolution 
with Vtg is influenced by both Ap and Ctg- Conse¬ 
quently, a simultaneous fit allows for a reliable estimation 
of Ap (through the doping dependence of ATg(Vtg)), 
and, in turn of Ctg (through the slope of Awg(Vtg) 
curve, having Ap constrained by ATg(Vtg))- 

Figure [5] shows the evolution of Ep as a function of 
Awg for five different graphene FETs (denoted sam¬ 
ple 1 to 5) in which Ctg and Ar have been previously 
determined by the simultaneous fit of Awg(Vtg) and 
ATg (Vtg)- For these five samples, we found an av¬ 
erage of (Ctg) = (4-5 ± 1.5) /iF cm” 2 and (Ar) = 
(0.032 ± 0.004) (see also Table Q. This translates into 
an average relative G-mode FWHM (see Eq. |[5|) of 
ATg = 12.6 ± 1.6 cm” 1 (at T = 0, Ep = 0 and SEp = 0) 
that is consistent with the value of Eg ~ 15 cm” 1 
recorded on quasi-undoped suspended graphene at low 
temperature. 77 Remarkably, and in spite of the different 
values of Ctg, the data for these five devices shown in 
Fig. [5] collapse onto a same curve. In practice, this very 
reproducible behavior can be used to evaluate Ep know¬ 
ing Awg, which is of broad interest in graphene science. 
For this purpose, we consider the asymptotic behavior of 
Aujg(Ep). When |Cp| Hcoq/2, Eq. 0 becomes 
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Assuming that the adiabatic contribution A ioq is neg¬ 
ligible compared to Acc2 A , Awg should be linear with 
\E f \. Indeed, in Fig. [5] for the five different samples, 
wg(^f) clearly scales linearly for \E F \ > 100 meV. The 
slightly different slopes observed for electron and hole 
doping arise from the opposite sign of the adiabatic cor¬ 
rections. 

For |Sp| > 100 meV, we find 


Ep > +100 meV, E F = +21Aw G + 75, (9) 

Ep < -100 meV, E F = -18Aw G - 83, (10) 


where Ep is expressed in meV and Aw G in cm -1 . How¬ 
ever, it should be noted that this linear scaling only holds 
for \E F \ < 500 — 600 meV. In fact, for higher E F \. 
Aw g no longer scales linearly with |.Ep| since Acan 
no longer be neglected compared to Awq A . Moreover, 
Eqs. (|9| and (10) can be applied provided the shift in 
Aw g is exclusively due to doping, i.e., other extrinsic 
factors, such as mechanical strain do not contribute. If 
this is not the case, one has to separate the various con¬ 
tributions, using, e.g. the method described in Ref. [78] 
with the results of Sec. IVII Bl 


B. Doping-dependence of the 2D-mode feature 

Let us briefly comment on the 2D-mode feature. Fig¬ 
ure [6] shows the evolution of the frequency w 2 d(Ef) and 
FWHM r 2 D (L?f) °f the 2D-mode feature with Ep for 
sample 2 (2D mode spectra are also shown for sample 
1 in Fig. [2]). In supported graphene, the 2D-mode fea¬ 
ture typically exhibit a quasi-symmetric lineshape that 
can b e phen omenologically fit to a modified Lorentzian 
profile We find that F 2 d(7?f) does not vary sig¬ 
nificantly with the gate bias, while 102 d varies little at 
moderate doping (|J?f| < 200 meV), but tends to stiffen 
(soften) significantly for stronger hole (electron) doping. 
The observed evolution of w 2 d outlined in Fig. [6] (see also 
Figs. [ 2 ] and |9]) can be qualitatively understood as the sum 
of a dominant adiabatic contribution and a weaker non- 
adiabatic contribution. The latter is reduced as com¬ 
pared to the case of the G-rnode feature, likely because 
the 2D-mode feature involves phonons that are signifi¬ 
cantly away from the edges of the Brillouin zone.® 


n (x 10 12 cm' 2 ) 
0.0 2.4 



FIG. 6. Doping dependence of the frequency cu 2D ( re d squares, 
left axis) and FWHM F 2D (blue circle, right axis) of the 2D 
mode-feature. The measurements are performed on sample 2, 
before the creation of defects. 


E F j 15 l 41 l 44 l whereas / G does not, as long as \Ep\<Jud]j2, 
where hio F is the energy of the incident laser.®® 79] j n 
Fig. 0 we consider the ratio / 2 d/ 7 G , which is maximum 
for E F = 0 and decreases almost symmetrically for in¬ 
creasing \E F \. Following Ref. U5l and 1551 the integrated 
intensity of the 2D-mode feature writes 


h d oc 


7k 


7e— P h + 7D + 7e. 


( 11 ) 


where 7 e - P h + 7d + 7ee is the total electron scattering 
rate, with 7 e - P h the electron-phonon scattering rate, 7 d 
the electron-defect scattering rate, and 7 ee the electron- 
electron scattering rate. The electron-phonon scattering 
rate can be approximated as 7 e - P h = 7k + 7r, where 7k 
and 7r are the scattering rates for zone-edge and zone- 
center optical phonons, respectively. Note that Eq. (11) 


is obtained under the assumption of a fully resonant 
process (see Fig. [2j[f)), and that trigonal warping ef¬ 
fects leading to momentum-dependent scattering rates 
are neglected! 15 l 59 l 6 ° l 8 Sl While 7 d and 7 G _ p h do not de¬ 
pend on Ep, 7 ee has been predicted to scale linearly with 
\E F \. For |Ek| < hidp/2, Basko et al. calculated^ 


C. Electron-electron and electron-phonon 
scattering 




7e- P h + 7D 


(7e- P h + 7D +0.06|-E f |), 


( 12 ) 


Another useful quantity is the integrated intensity of a 
Raman feature (denoted lx), which represents the total 
probability of the Raman scattering process. The inte¬ 
grated intensity of the 2D-mode feature (/ 2 d) depends on 


where 


corresponds to the value at Ep = 0. 


In this section, we are considering pristine graphene, 
in which <C 7 e - P h- As illustrated by the dashed line 
in Fig. [7] our experimental data agree well with a fit 
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FIG. 7. Left axis: doping dependence of the ratio between 
the integrated intensity of the 2D-mode feature and that of 
the G-mode feature. Right axis: doping dependence of the 
square root of the ratio between the integrated intensity of 
the G mode and that of the 2D mode. The dashed and solid 
lines are fits based on Eq. (121, without and with Fermi en¬ 
ergy fluctuations, respectively. A value of 7 e - P h = 51 meV is 
deduced from the fit. The measurements are performed on 
sample 2, before the creation of defects. 


based on Eq. (12 1 for |£Jp| > 100 meV. However, we 
observe a deviation from Eq. (12) near the CNP, likely 
due to Fermi energy fluctuations. As in Sec. V A| we 
therefore fit the experimental data with the Gaussian 
convolution of Eq. (12), resulting in the solid line in 
Fig-0 The agreement between theory and experiment 
is very good and more compelling than in the seminal 

study in Ref. dH The fitting parameters are 7 P =3.6, 

G o 

7e-ph = 51 meV and 6E F = 110 meV. We repeated this 
analysis on three pristine samples and found average val¬ 
ues of { 7 e_ ph } = (47 ± 7) meV, (SE F ) = (120 ± 10) meV 


and 


12D 

Ig 


= 4.2 


I). Note that 
lese three de- 


± 0.6 (see Table 

the dispersion of the measurements on t 
vices is very similar to the dispersion observed when 
measuring on several spots on the same sample. The 
value of 7 e_ p h is in good agreement with the estimate in 
Ref. HU The Fermi energy fluctuation SEp obtained here 
is more realistic than the lower values estimated from 
the simul taneous fit of Ao>g(Ftg) and AFg(Vtg) (see 
Sec. VA). It corresponds to a charge inhomogeneity of 
Sn < 10 l2 cm' 2 , in line with previous scanning tunneling 
microscopy measurements EALZAl 

Interestingly, in Ref. |T5] the authors claim that the 

intrinsic value of for undoped graphene is in 

G o 

the range 12-17 (using a 514.5 nm excitation wave¬ 
length). However, this estimation is based on Raman 


measurements on quasi-undoped suspended graphene 19 
and does not take into account the effect of optical inter¬ 
ferences, which occur in graphene-based multilayer struc¬ 
tures and may critically affect the intensity of the Ra¬ 
man features! 81 * 82 ! From the data in Ref. [H2J an intrin¬ 


sic value corrected from interference effects of 


1 2D 

Ig 


intr 


5 ± 0.3 can be estimated for freely suspended, undoped 
graphene, using a 532 nm excitation wavelength, as in 
the present study. Considering the distinct Raman en¬ 
hancement factors for the G- and 2D-mode features in 
the PEO/graphene/Si 02 (285 nm)/Si multilayer struc¬ 


ture, our average value of 


^2D 

Ig 


= 4.2 ± 0.6 translates 


into an average intrinsic value of ^ y = 4.9 ± 0.7 

(see Table [I]), which is in excellent agreement with our 
estimate on suspended graphene. 

As outlined in Ref. nasi, and (S 3 , the scattering rate 
7e_ p h is linked to the dimensionless electron-phonon cou¬ 
pling constants Ar and \k through 


7e-ph = 7K+7r 


Ar { TlU L 

~T l ~Y~ 


huiK 



fiwL 

~Y~ 



(13) 

where Hljk ~ 1210 cm 1 = 150 meV is the in-plane 
transverse optical (TO) phonon energy at the K (K’) 
point, Hlot '■= Eo g ~ 1580 cm ' 1 = 196 meV is the in¬ 
plane optical phonon energy at T (i.e., the G-mode fre¬ 
quency) and fno l = 2.33 eV is the laser photon energy. 

For sample 2 (see Fig. [7] and Table [T]), a value of 
Ar = 0.034 is deduced from the simultaneous fits of 
Aojg and ATg (see Sec. |V Aj ). Then, using Eq. (13), we 
can estimat(pAl / \ K = 0.17. Overall, for the three pris¬ 
tine samples studied here, we obtained average values 
of (A r ) = 0.031 ± 0.004, (A K > = 0.15 ± 0.03 and 

= 5.1 ± 1.2 (see Table|lj). 

To close this section, we compare the average ra¬ 
tio deduced from our doping-dependent Raman 


study to a direct estimate derived from the measured 
ratio of the integrated intensities of the intravalley (2D’ 
mode) and intervalley (2D mode) resonant two-phonon 
features F 11 This ratio is expected to be independent of 

Ep and writeipEIlS jm. = 2 (^ j . In our experimen¬ 
tal conditions, we obtain =40 ± 2. Thus, by 
considering one more time the different Raman enhance¬ 
ment factors for the 2D- and 2D’-mode features in the 
PEO/graphene/Si 02 (285 nm)/Si multilayer system, we 
deduce ^ ss 4.2. This value is consistent with the anal¬ 
ysis outlined above. 


VI. DEFECTIVE GRAPHENE 
A. Creation of defects 

As mentioned in Sec. [TTJ when an electrochemically 
gated graphene FET is subjected to a sufficiently high 
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Sample 

J 2P 

Ig 

intr 

(fe)„ 

n D (xl0 12 cm 2 ) 

7e- P h + 7 d (meV) 

A r 

Ak 

1 (without defects) 

5.6 

< 0.05 

- 

50 

0.027 

0.17 

1 (with defects) 

4.6 

1.7 

0.9 

57 

0.031 

< 0.20 

2 (without defects) 

4.2 

< 0.05 

- 

51 

0.034 

0.17 

2 (with defects) 

3.3 

1.3 

0.7 

69 

0.031 

< 0.24 

3 

5.0 

< 0.05 

- 

39 

0.031 

0.12 

4 

4.4 

2.6 

1.4 

53 

0.037 

< 0.18 

5 

4.3 

1.4 

0.7 

72 

0.031 

< 0.25 


TABLE I. Intrinsic integrated intensity ratio (corrected from interference effectJ^) 1 at Ep = 0, measured integrated 


intensity ratio 


i 

!gJ 




at Ep = 0, estimated defect concentration «d, sum of the electron-phonon and electron-defect scattering 


rates, and dimensionless electron-phonon coupling constants at T and near K (K’), for five different electrochemically gated 
graphene transistors. The measurements be fore and after the creation of defects have been done (i) at the same spot in sample 
1 and correspond to the data shown in Fig. 10 and (ii) at two different spots in sample 2. 


gate bias, electrochemical reactions may occui^^mi an( f 
create defects in the graphene channel. In our devices, 
a reaction systematically occurs at negative gate biases 
(Vt G k -1 V to -2 V). The threshold voltage depends 
on the sample and on the gate capacitance. Electrochem¬ 
ical reactions result in an increase of the gate leak cur¬ 
rent above 1 nA, and in the emergence of defect-induced 
features in the Raman spectrum. Figure [8] shows two Ra¬ 
man spectra recorded at Vtg = 0 V on sample 1, before 
applying any gate voltage and after an electrochemical 
reaction has taken place. We clearly see that (i) the G- 
and 2D-mode features do not shift, and (ii) prominent D- 
and D’-mode features develop. These two Raman modes 
are known to be forbidden by sym metry an d ca n only be 
observed in the presence of defects. ®^ 1 ' | l 62 l 85 l Hni 


B. Doping dependence of the Raman features in 
defective graphene 


Figure [9])a) shows wg(A’f) and ATg(-Ff) in defective 
graphene. By comparing this figure with Fig. [4j we con¬ 
clude that the doping dependence of the G-mode feature 
is not affected by the presence of defects. Both Aluq(E f ) 
and ATg(-Ff) are well fit to the theoretical model of 
Sec. EB The frequencies W 2 P and wp are also shown 
as a function of E F in Fig. [9](b). The D- and 2D-mode 
features are fit to a modified Lorentzian profileP®”^ We 
note that both frequencies follow identical trends. More 
precisely W 2 P ~ 2wp — 6 in cm -1 . The factor 2 is ex¬ 
pected, since the 2D mode is the two-phonon overtone of 
the D-mode (in the D-mode process, one inelastic scat¬ 
tering by a near zone-edge TO phonon is replaced by an 
elastic scattering by a def’edP ^ 60 * 86 ). This small differ¬ 
ence between 2and W 2 P is consistently observed in all 
the studied samples. It coul d be d ue to slight differences 
in the resonance conditions! 11 ® 87 



Raman shift (cm 1 ) 

FIG. 8. Raman spectra at Vtg = 0 recorded at the same 
location on sample 1, before any gate bias has been applied 
(black line) and after an electrochemical reaction has occurred 
(red line). The spectra are vertically offset for clarity. The 
dotted lines correspond to the baseline. 


C. Electron-defect scattering 


In Sec. VC using Eq. (12), we have shown that it 


is possible to deduce the phonon scattering rate 7 e - P h 
in pristine graphene from the study of the integrated 
intensities of the G- and 2D-mode features. In defec¬ 
tive graphene, a similar analysis can be performed pro¬ 
vided that a finite electron-defect scattering rate qp, pro¬ 
portional to the d efect concentration np, is taken into 
consideration .EDMI 

"a), we plot J as a function of E F . The 


In Fig. 

data is extracted from another series of measurements in 
sample 1, at a same spot, before and after the creation of 
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FIG. 9. Doping dependence of the Raman features in sample 
1, after the creation of defects, (a) Doping dependence of the 
frequency lo g (red squares, left axis) and relative FWHM Apjg 
( blue circles, right axis) of the G-mode feature. The blue and 
red solid lines are fits based on Eq. § and §, as in Fig. PH 
respectively. The fitting parameters are Ctg = 3.4 /rF cm - , 
Ar = 0.035 and SEp = 30 meV (b) Doping dependence of 
the frequencies of the 2D- (red squares) and D-mode (blue 
open circles) features (denoted W 2 D and ojd, respectively) in 
defective graphene. 


defects. As expected, we observe that the two datasets 
are well fitted by Eq. (12). The fitting parameters are 


7e- P h = 50 meV, ^ 


no JJ 


= 4.8, SE f = 110 meV for 


pristine graphene, and 7 e - P h+7D = 57 meV, 

SEp = 130 meV, for defective graphene, respectively. 


= 4.0, 


Since 7 e - P h is not affected by the presence of defects, 
we can estimate that 7 d ~ 7 meV for sample 1. Thus, 


although the D- and 2D-mode features have similar inte¬ 
grated intensities (Ip// 2 p ~ 0.5, as shown in Fig.[Io](b)), 
7d 7e- P h> as predict ed in t he low-defect concentration 


regime^ (see also Sec. VID). Another way to determine 


in the presence 


7 d is to compare the quantity ■J j ^ 

and in the absence of defects. According to Eq. (11) and 
(12), one obtains 


= 1 + 


The 


results 

D 


of Fig. 10 a) 


> 


7d 

7e-ph 

indeed show 
in agreement 


(14) 

that 

with 


Eq. (|14j). From the fitting parameters, we estimate 
that 7 d ~ 5 meV, in good agreement with the other 
estimate obtained above. 

To conclude this subsection, we focus on the de¬ 
pendence of Id on Ep. In practice, Ip is r outinely 
used to estimate a defect concentration. ilEHstKIl] jj ow _ 
ever, although it is not fully resonant, the D mode 
ma y involve o ne resonant electron-phonon scattering pro- 
ce S jllM H 62 l 92 l (see Fig. [2j)g)) - In other words, similarly to 
/ 20 , Id is also expected to decrease with increasing Ep \. 


In Fig. 10'b), we show jff and as a function of 

c) displays ^ and as a f unC - 

y, Jp and I 2 D show a very similar dop- 


Ep, while Fig. 


10 


tion of Ep. Clear, 
ing dependence (see also Ref. [HD ■ More quantitatively, 

a phenomenological fit of y^(Ep) using Eq. (12) (ap¬ 
plied to Id instead of / 20 ) agrees well with our measure¬ 
ments (see Fig. 10 'c)) and yields 7 e - P h + 7 d =52 meV, 


In 

Ig 


= 2.2, and SEp = 130 meV. The value of 7 e - P h+7D 


is very close to that obtained by fitting w > 

We note that Ref. m report a value of 7 e - P h + 7d 
70 meV similar to ours (see Table [I]), for defective 
graphene samples with slightly larger, yet similar val¬ 
ues of 4°- . However, a larger value of qp ~ 40 meV 
G 0 

is estimated, using the value of 7 e - P h ~ 30 meV ex¬ 
tracted from the measurements in Ref. EH. Our work 
provides an estimate of yp from a series of measure¬ 
ments performed on a same sample and suggests that 

7d 7e—phi even for > 1. Consequently, in 

lG 0 

defective graphene, we consider that the slope of the 
yjjf^(Ep) curve provides a fair estimate of 7 e - p hi from 
which we extract Ak, knowing Ap. Albeit the existence 
of a finite yp presumably leads to a slight overestimation 
of Ak, we do not observe a large difference between the 
values measured on defective and on pristine graphene 
(see Table [i]) . More quantitatively, by averaging on four 
defective graphene samples with similar defect concen¬ 
trations, we obtain (7 e - P h + 7d) = (63 ± 9) meV, 
a value that is indeed slightly larger than the average 
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FIG. 10. (a) Integrated intensity ratio y as a function of Ef in sample 1 before (red open squares) and after (red circles) the 

creation of defects, (b) Integrated intensity ratios ^ (red circles), (green triangles) and ^ (blue squares), as a function of 

Ef in sample 1, after the creation of defects, (c) \J (red circles, left axis) and (blue squares, right axis) as a function 

of Ef, in sample 1, after the creation of defects. As described in Sec. VID a defect concentration tid 
estimated from these measurements. 


(7e_ p h) = (47 ± 7) meV obtai ned on three pristine sam¬ 
ples (see Table |T] and Sec. VC). 


10 14 

"D = -J2- 


1.8 x 10 10 (hu h ) 4 


(15) 


D. Defect concentration 

In principle, the concentration «d of defects in a 
graphene sample can be deduced from the analysis of 
the defect-induced Raman modes, such as the (interval¬ 
ley) D mode or the (intravalley) D’ mode. The study 
of the defect-induced Raman modes has far reaching 
consequences for sample characterization and can also 
be a very useful tool to monitor chemical reactions on 
graphene. Following the seminal work by Tuinstra and 
Koenigp^ several groups have proposed analytical ex¬ 
pressions to connect Id and tid in various graphitic ma¬ 
terials, from weakly defective graphene layers to amor¬ 
phous carbonfS^HESI 

The defective graphene samples studied here exhibit 


where «d is the concentration of defects in cm -2 , Ld is 
in nm, Hoj l is the laser photon energy in eV and ^ 

is taken at Ef ~ 0, still with Fermi energy fluctuations. 
According to Ref. 1111 the scaling introduced in Eq. (15) 
is independent of the type of defect. 

In Fig. till we plot jP (Ef), normalized by its value at 
Ef ~ 0 for five different sets of measurements (including 
some at different locations on the same sample), with 
different defect concentrations. We observe that all the 
data collapse onto the same curve. For \E F \ < 100 meV, 

7P ss (?g) anc ^ this ratio decreases by less than 20 % 

for |Ep| < 200 meV. Thus, since unintentional doping in 
graphene samples typically leads to \Ef\ < 200 meV, the 


an integrated intensity ratio r ° < 2.6 near the CNP experimentally measured can be used together with 


(see Fig. H and [TTJ and Tabic Their Raman features 
show only a slight spectral broadening (by a few cm -1 ) 
compared to the pristine case (see Fig. [8]). More pre¬ 
cisely, we find, for five series of gate-dependent measure¬ 
ments on various defective regions, that To ~ 10 cm -1 
(as compared to To ~ 5 cm -1 for pristine graphene), 
T2D ~ 30 cm -1 and Td ~ 20 cm -1 . Following the three 
stage classification of Ref. [55] and related works p4S2J EH 
such samples can be described as stage 1, i.e., still in 
the weakly defective regime. Let us assume point de¬ 
fects, separated by an average distance Ld 10 nm. In 
this regime, Eq. (9) of Ref. [ST] and the results of Ref. 021 
provide the relational 


Eq. (15) for an estimation of no in weakly doped samples. 


Using Eq. (15), the electrochemically-induced defect 


concentrations d edu ced for our measurements (see Ta¬ 
ble [I] and Fig. 11) range from 2.7 x 10 11 cm -2 to 


1.4 x 10 12 cm -2 . This translates into Ld ranging from 
19.5 nm down to 8.5 nm. The latter value is at the limit of 
the weakly defective regime, which assumes Ld > 10 nm. 

Overall, the results shown in Figs. m demonstrate 
that for defect concentrations below approximately 2 x 
10 12 cm 2 , the electron-defect scattering rate remains 
much smaller than the electron-phonon scattering rate, 
and that the doping dependence of the G- and 2D-mode 
features is essentially the same as in pristine graphene. 
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FIG. 11. Doping dependence of the integrated intensity ratio 
for five sets of measurements on defective graphene sam¬ 
ples. Each dataset is normalized by the value (tJt) mea¬ 
sured near the charge neutrality point. The measured (j 
are indicated in the legend. 


These results contrast with the fact that even for rela¬ 
tively low np> in the range 10 11 — 10 12 cm -2 , the inte¬ 
grated intensity of the D-mode feature is smaller, yet 
on the same order of magnitude as that of the 2D- 
mode feature, in keeping with recent experimentaP 6 87 1 
and theoretical resultsPS This calls for further investi¬ 
gations of the integrated intensity of the one-phonon, 
defect-induced Raman features relative to that of their 
symmetry-allowed overtones. 

VII. CORRELATIONS 

In the previous sections, we have successfully com¬ 
pared our measurements to theoretical calculations and, 
in particular estimated the electron-phonon coupling con¬ 
stants. In this section, we present correlations between 
the frequencies and linewidths of the main Raman fea¬ 
tures in doped graphene, with the aim to extract uni¬ 
versal behaviors that could be useful for sample char¬ 
acterization. Based on the conclusions of Sec. m the 
correlations discussed in the following will also hold in 
weakly defective graphene. 

A. G-mode frequency and linewidth 

Figure [12| shows ATg as a function of A ujq for the 
five different samples already shown in Fig. [5j We ob¬ 
serve a universal behavior and the experimental data are 



FIG. 12. Correlation between the relative FWHM AFg 
and relative frequency Aojg of the G-mode feature in doped 
graphene, for the five samples introduced in Fig. [5] The solid 
and dashed lines correspond to theoretical calculations (based 
on Eqs. © and ©), for electron and hole doping, respec¬ 
tively. 


in good agreement with the theoretical calculations, al¬ 
though the very slight difference expected for electron 
and hole doping (due to AWq, see Eq. <©) is not re¬ 
solved experimentally, likely due to Fermi energy fluc¬ 
tuations. We also note that in the high-doping regime, 
ATg tends to increase somewhat. This increase, also ob¬ 
served by others is presumably due to the increasing 
inhomogeneity of the charge distribution at high top-gate 
biases. The correlation displayed in Fig. [l2|may also be 
used to estimate Ep, especially in the low doping regime 
(|£ F |<fit4). 

B. G- and 2D-mode frequencies 

Figure [13] represents the evolution of W2D as a function 
of ug for the same five samples. A clear correlation is 
observed between these two quantities. For hole doping, 
the correlation is quasi-linear in the range of Ep stud¬ 
ied here (—500 meV < Ep < 0). In contrast, for elec¬ 
tron doping, a quasi-linear scaling, again with a (much 
smaller) positive slope is also observed at low doping 
(0 < Ep < 250 meV), until W2D levels off and ultimately 
decreases, leading to a non-linear scaling. This behavior 
was observed on every sample either for electrolyte-gated 
or conventional back-gated FETs and has been also ob¬ 
served in chemically doped graphene. 78 

From the slopes extracted on approximately 

thirty samples (see Fig. M b ) -(c)), we find an average 
of (0.55 ± 0.2) for hole doping and of (0.2 ± 0.13) for 
electron doping, respectively. The former value agrees 
well with the slope of (0.70 ± 0.05) extracted numer¬ 
ically by Lee et all 78 from the data in Ref. [JI] and [H] 
From our statistical study, we note that the correlation 
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FIG. 13. (a) Correlation between the frequencies of the 2D- 
and G-mode features (relative to their values near the charge 
neutrality point) in doped graphene, for the five samples in¬ 
troduced in Fig. [5] The continuous and short-dashed lines are 
global fits performed on the linear portions of the hole dop¬ 
ing and electron doping branches, respectively. The dashed 
line corres ponds to the evolution of a> 2 D versus cjg under 
pure strain!- 78 ^! Statistical distribution of the measured slopes 
under hole and electron doping are shown in (b) and (c), 
respectively. 


between W2D and wq is more dispersed than the correla¬ 
tion between A1?g and Awe- This is chiefly due to the 
dependence of uj 2D on Ep, which is not as universal as 
that of U(j. In addition, it is rather challenging to extract 
a well-defined correlation for electron doping due to the 
small variations of W2D at moderate doping. 

Noteworthy, estimations of Ep based on the frequency 
and/or linewidth of the Raman features may only be re¬ 
liable if graphene is not subjected to significant strains. 
Indeed, Tq is marginally affected by isotropic strains be¬ 
low l°/cP However, the G-mode feature may broaden 
and ultimately split into two su b-feat ures in the pres¬ 
ence of larger anisotropic strains p95| j n addition, the 
Raman features soften (stiffen) under tensile (compres¬ 
sive) strain. A linear correlation between W2D and uq 
has been measured in strained graphene ! 78 * 82 * 96 * ^ Since 
the measured slopes (~ 2.2, for undoped, strained 
graphene®) are appreciably larger than the slopes mea¬ 
sured in doped graphene (presumably under a small but 
constant built-in strain), Lee et al. have proposed to use 



FIG. 14. Charge carrier density n as a function of the shift 
[(ug — c4) 2 + (^2D — W2 D ) 2 ] 1 ^ 2 from the reference point cor¬ 
responding to undoped graphene. Data are shown for the five 
samples introduced in Fig. [5] The continuous and dashed 
lines are linear global fits performed on the electron and hole 
branches, respectively. 


the correlation between W2D and wg as a robust tool to 
optically separate strain from charge dopingP 

Following Ref. 178] we may then define three vectors 
corresponding to the slopes under strain, hole and 

electron doping, respectively (see Fig. 131. To further 
deduce absolute levels of strain and/or doping, one also 
has to know the 2D- and G-mode frequencies that corre¬ 
sponds to an undoped and unstrained graphene sample. 
For clarity, in Fig. [l3j the 2D- and G-mode frequencies 
are shown relative to the measurements at Ep ss 0. These 
origin points, denoted (wq, w 9 d ) might differ from the 
reference point corresponding to undoped and unstrained 
graphene, since an undetermined amount of native strain 
may be present and induce a shift along the strain vector. 


The data in Fig. fTsl allows an estimation of the coef- 
ficient which connects [(wq -Wq) 2 + (w 2 d — w 9 D ) 2 ] j 
the measured distance from the zero doping point, to a 
given doping level (see Fig. [l4]). We chose to consider n 
instead of Ep because it is a more relevant quantity as 
far as graphene characterization is concerned. Although 
the curves displayed in Fig. [14] are not expected to ex¬ 
hibit a linear scaling (as opposed to the data shown in 
Fig-©, we observe a quasi-linear scaling for sufficiently 
small doping (|n| <7x 10 12 cm -2 ). We therefore fit the 
linear part for both electron- and hole-doping with a line 
intercepting the zero doping point. We find slopes of 
4.4 x 10 11 cm -1 for electrons and —3.6 x 10 11 cm -1 for 
holes, respectively. 


Finally, considering the Griineisen parame¬ 
ters of 1.8 and 2.4 for the Raman G- and 2D- 
modes under isotropic strain,^! we can estimate 
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a slope of 7.1 x 10 3 % strain/cm 1 to connect 

[(wq — w g) 2 + ( W 2 D — w °d) 2 ] ~ t° an applied isotropic 

strain. In practice, the strain field may be anisotropic, 
depending on the sample and on the experimental condi¬ 
tions, leading to a different slope. These coefficients may 
be used for a reliable estimation of doping and strain in 
graphene samples and devices. 

VIII. CONCLUSION 

We have presented a robust method, based on Raman 
scattering spectroscopy, to accurately determine the ge¬ 
ometrical capacitance, and hence, the Fermi energy in 
electrochemically-gated graphene field-effect transistors 
with a spatial resolution down to approximately 1 /xm. 
Such a calibration allows for quantitative analysis of the 
doping dependence of the frequency, linewidth and in¬ 
tegrated intensity of the main Raman features. The 
anomalous doping dependence of the G-mode phonons is 
well captured by theoretical models over a broad range 
of Fermi energies above or below the Dirac point, and 
provides an experimental measurement of the electron- 
phonon coupling constant at the T point of the Brillouin 
zone. We have then exploited the peculiar doping de¬ 
pendence of the integrated intensity of the multiphonon 
resonant Raman features, in particular the resonant 2D- 
mode feature, to estimate the electron-phonon coupling 
constant at the edges (K, K’) of the Brillouin zone. Fi¬ 
nally, from the doping dependence of the integrated in¬ 
tensity of the defect-induced D-mode feature, we can es¬ 
timate the electron-defect scattering rate in stage 1 de¬ 
fective graphene samples. 

Our study provides useful guidelines for the character¬ 
ization of graphene samples. We have, in particular, con¬ 
sidered the correlation between the frequency and width 
of the G-mode feature, as well as between the frequencies 
of the 2D- and G-mode features. These correlations re¬ 
veal universal behaviors that can therefore be applied to 


evaluate doping in a variety of experimental situations. 
We have also demonstrated that defects can be efficiently 
created in-situ in electrochemically gated graphene field 
effect transistors. The integrated intensity of the D-mode 
feature decreases monotonically with increasing doping, 
and follows the same scaling as that of its two-phonon 
overtone. However, due to Fermi energy fluctuations, 
the D-mode intensity is nearly constant for Fermi energy 
shifts below 200 meV relative to the Dirac point, which 
is of practical interest for the determination of the defect 
concentration. 

In the present work, we conservatively estimate that 
Fermi energies as high as ~ 700 meV above the Dirac 
point can be achieved in ambient conditions, without 
damaging graphene. This naturally opens exciting per¬ 
spectives for optoelectronics. Nevertheless, a well con¬ 
trolled Raman scattering study of the crossover be¬ 
tween the intermediate doping regime achieved here 
('n ss 10 13 cm -2 ) and the very high doping regime 
(n > 10 14 cm -2 ) remains elusive. Finally, electrochemi¬ 
cal gating is a promising strategy to investigate electron- 
phonon coupling in other two-dimensional materials, in¬ 
cluding transition metal dichalcogenidesP^ 
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